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Abstract 

The first observation of the decay B® — > D°K~ir + is reported. The analysis is 
based on a data sample, corresponding to an integrated luminosity of 1.0 fb" 1 of 
pp collisions, collected with the LHCb detector. The branching fraction relative to 
that of the topologically similar decay B° — > D tt + it~ is measured to be 

B (B° s -> D°K-7r+) , , 

-A^ { = 1.18 ± 0.05 stat. ± 0.12 syst. . 

In addition, the relative branching fraction of the decay B° — >• D°K + tt~ is measured 
to be 

B (B° -> D°K+tv-) , . , . 

-±— , ; = 0.106 ± 0.007 (stat.) ± 0.008 (syst.) . 

B (B U -> Z? U 7T + 7T~) 
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1 Introduction 



The precise measurement of the angle 7 of the CKM Unitarity Triangle [HE] is one of the 
primary objectives in contemporary flavour physics. Measurements from the experiments 
BaBar, Belle and LHCb are based mainly on studies of B + — > DK + decays, where the 
notation D implies that the neutral D meson is an admixture of D° and D° states. Each 
experiment currently gives constraints on 7 with a precision of ~ 15° [3H5]. Significant 
reduction of this uncertainty is well motivated and the use of additional channels to further 
improve the precision is of great interest. 

The decay B° — > DK + w~, including the resonant contribution from B° — > DK*°, is 
one of the modes with the potential to make significant impact on the overall determi- 
nation of 7 [6] . A first measurement of CP observables in B° — > DK*° decays has been 
reported by LHCb [7]. This decay is particularly sensitive to 7 owing to the interference 
of b — > cus and b — > ucs amplitudes, which for this decay are of similar magnitude. It has 
been noted that an amplitude analysis of B° — > DK + ir~ decays can further improve the 
sensitivity and also resolve the ambiguities in the result [HIE]- 

The decays B° — > D°K + ir~ and B° s — > D°K~ir + can be mediated by the decay dia- 
grams shown in Fig. [TJ Both B° and B° s decays are flavour-specific, with the charge of the 
kaon identifying the flavour of the decaying B meson, though the charges are opposite in 
the two cases. In addition to these colour-allowed tree-level diagrams, colour-suppressed 
tree-level diagrams contribute to B? g s — > D°Kir decays (Kir denotes the sum over both 
charge combinations). Both colour-allowed and colour-suppressed diagrams contribute to 
the CKM-suppressed B9 \ — > D°Kii modes. 

A first study of the decay B° — > D°K + ir~ has been performed by BaBar [TU], giving 
a branching fraction measurement B — >■ D°K + n~) = (88 ± 15 ± 9) x 10~ 6 , where 
the contribution from the B° — >■ D*~K + decay is excluded. There is no previous branch- 
ing fraction measurement for the inclusive three-body process B® — > D°K~ti + , although 
that of the resonant contribution D°K*° has been measured by LHCb [TT]. Since the 
B® — > D°K~ii + and the related B® — > D*°K~7i + decays form potentially serious back- 
grounds to the -B — > D K + ir~ channel, measurements of their properties will be neces- 
sary to reduce systematic uncertainties in the determination of 7. 

In this paper the results of a study of neutral B meson decays to D°Kii, includ- 




Figure 1: Decay diagrams for (a) favoured B° — > D°K + tt decays and (b) favoured 
B° -)• D°K~ir + decays. 
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ing inspections of their Dalitz plot distributions, are presented. The D°K + ir and 
D°K~tt + final states are combined, and the inclusion of charge conjugate processes is 
implied throughout the paper. In order to reduce systematic uncertainties in the mea- 
surements, the topologically similar decay D°n + ix~, which has been studied in detail 
previously P^IES], * s used as a normalisation channel. In this paper, Dim denotes the 
D°7r + 7r~ final state and DK7r denotes the sum over the D°K + ir~ and D°K~tt + final 
states. The neutral D meson is reconstructed using the D° — > K + n~ final state; therefore 
the signal yields measured include small contributions from D° — >■ K + n~ decays, but such 
contributions are expected to be small and are neglected hereafter. The analysis uses a 
data sample, corresponding to an integrated luminosity of 1.0 fb" 1 of pp collisions at a 
centre-of-mass energy of 7TeV, collected with the LHCb detector during 2011. 

2 Detector, trigger and selection 

The LHCb detector [H] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detec- 
tor includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4Tm, and three stations of silicon- 
strip detectors and straw drift tubes placed downstream. The combined tracking system 
provides a momentum measurement with relative uncertainty that varies from 0.4% at 
5 GeV/c to 0.6 % at 100 GeV/c, and impact parameter (IP) resolution of 20 [im for tracks 
with high transverse momentum (pt)- Charged hadrons are identified using two ring- 
imaging Cherenkov (RICH) detectors [15]. Photon, electron and hadron candidates are 
identified by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a 
system composed of alternating layers of iron and multiwire proportional chambers. 

The LHCb trigger [16] consists of a hardware stage, based on information from the 
calorimeter and muon systems, followed by a software stage that applies a full event 
reconstruction. In this analysis, signal candidates are accepted if one of the final state 
particles created a cluster in the hadronic calorimeter with sufficient transverse energy 
to fire the hardware trigger. Events that are triggered at the hardware level by another 
particle in the event are also retained. 

The software trigger requires a two-, three- or four-track secondary vertex with a high 
sum of the transverse momentum, p^, of the tracks and a significant displacement from 
the primary pp interaction vertices (PVs). At least one track should have > 1.7GeV/c 
and impact parameter x 2 , Xip> with respect to the primary interaction greater than 16. 
The Xip is the difference between the \ 2 °f the PV reconstruction with and without the 
considered track. A multivariate algorithm [T7] is used for the identification of secondary 
vertices consistent with the decay of a b hadron. 

Candidates that satisfy the software trigger selection and are consistent with the decay 
chain B?, — > D°K ± ti^, D° — >■ K + tt~ are selected, with requirements similar to those in 
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the LHCb study of the decay B9-. — > D°K + K [IB]. The D° candidate invariant mass 
is required to satisfy 1844 < < 1884 MeV/c 2 . Tracks are required to be consistent 
with either the kaon or pion hypothesis, as appropriate, based on particle identification 
(PID) information primarily from the RICH detectors [15]. All other selection criteria 
were tuned on the D 7i + n~ channel. The large yield available for the B° — > D°7t + tt~ 
normalisation sample allows the selection to be based on data, though the efficiencies are 
determined using simulated events. In the simulation, pp collisions are generated using 
Pythia 6.4 [H] with a specific LHCb configuration [20] • Decays of hadronic particles are 
described by EvtGen [21] in which final state radiation is generated using Photos [22]- 
The interaction of the generated particles with the detector and its response are imple- 
mented using the Geant4 toolkit [23] as described in Ref. [24] . 

Loose selection requirements are applied to obtain a visible signal peak in the D 7t + tt~ 
normalisation channel. The selection includes criteria on the quality of the tracks forming 
the signal candidate, their p, p T and Xip- Requirements are also placed on the corre- 
sponding variables for candidate composite particles (D°, B? s -A together with restrictions 
on the consistency of the decay fit (Xvertox)' the flight distance significance (Xm g ht)' anc ^ 
the cosine of the angle between the momentum vector and the line joining the PV under 
consideration to the B^ vertex (cos#dir) [TT] . 

A boosted decision tree (BDT) [22] that identifies D° — > K + ir~ candidates is used to 
suppress backgrounds from fe-hadron decays to final states that do not contain charmed 
particles and backgrounds where the D° does not decay to the K + tt~ final state. This 
tl D° BDT" [2"6|l27] is trained using a large high-purity sample obtained from B + — > D°ti + 
decays. The BDT takes advantage of the kinematic similarity of all 6-hadron decays 
and avoids using any topological information from the B9-, decay Properties of the D° 
candidate and its daughter tracks, containing kinematic, track quality, vertex and PID 
information, are used to train the BDT. 

Further discrimination between signal and background categories is achieved by calcu- 
lating weights, using the sPlot technique [2H], for the remaining D°n + ir~ candidates. The 
weights are based on a simplified fit to the B candidate invariant mass distribution from 
the Dim data sample. The weights are used to train a neural network [21] to maximise 
the separation between the categories. A total of 10 variables are used in the network. 
They include the p?, yfp? Xvcrtex' Xm g ht anc ^ cos $dir of the B9-. candidate, the output 
of the D° BDT and the yfp °f the two pion tracks that originate from the B^ vertex. 
The pp asymmetry and track multiplicity in a cone with half- angle of 1.5 units in the 
plane of pseudorapidity and azimuthal angle (measured in radians) [30] around the B^ 
candidate flight direction are also used. The input quantities to the neural network only 
depend weakly on the kinematics of the B9 ^ decay. A requirement on the network output 
is imposed that reduces the combinatorial background by an order of magnitude while 
retaining about 70 % of the signal. 

To improve the B^ candidate invariant mass resolution, the four-momenta of the 

tracks from the D° candidate are adjusted [31] so that their combined invariant mass 
matches the world average value [32] • An additional B? s) mass constraint is applied in the 
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calculation of the Dalitz plot coordinates, m 2 (DK) and m 2 (Dir), which are used in the 
determination of event-by-event efficiencies. The coordinates are calculated twice: once 
each with a B° and a B° s mass constraint. A small fraction (~ 1 % within the fitted mass 
range) of candidates with invariant masses far from the B9 ^ peak fail one or both of these 
mass-constrained fits, and are removed from the analysis. 

To remove the large background from B° — > D*~ir + decays, candidates in both samples 
are rejected if the mass difference rriD^-mo (for either pion charge in the combinations 
D°7r + 7r~ and D°Ktt) lies within ±2.5MeV/c 2 of the nominal D*~-D° mass difference [32] . 
Candidates in the DKtt sample are also rejected if the mass difference m£)^-m£) calculated 
under the pion mass hypothesis satisfies the same criterion. A potential background 
contribution from B® — > D^K^ decays is removed by requiring that the pion from the D° 
candidate together with the kaon and the pion do not form an invariant mass in the range 
1850-1885 MeV/c 2 . Further DKir candidates are rejected by requiring that the kaon from 
the D° candidate together with the kaon and the pion do not form an invariant mass in the 
range 1955-1975 MeV/c 2 , which removes potential background from B° s — > D^w decays. 
A muon veto is applied to all four final state tracks to remove potential background from 
B® s j — > J/ip K*° decays and D° candidates are required to travel at least 1 mm from the 
B® s j decay vertex to remove charmless backgrounds that survive the D° BDT requirement. 

Candidates are retained for further analysis if they have an invariant mass in the 
range 5150-5600 MeV/c 2 for Dtxtx or 5200-5600 MeV/c 2 for DKir. After all selection re- 
quirements are applied, fewer than 1 % of events with at least one candidate also contain a 
second candidate. Such multiple candidates are retained and treated in the same manner 
as other candidates; the associated systematic uncertainty is negligible. 

3 Determination of signal yields 

The signal yields are obtained from unbinned maximum likelihood fits to the invariant 
mass distributions. In addition to signal contributions and combinatorial background, 
candidates may be formed from misidentified or partially reconstructed 6-hadron decays. 
Contributions from partially reconstructed decays are reduced by the lower bounds on 
the invariant mass regions used in the fits. Sources of misidentified backgrounds are 
investigated using simulation. Most potential sources are found to have broad invariant 
mass distributions, and are absorbed in the combinatorial background shapes used in the 
fits described below. Backgrounds from A b — > D°pir + [53] and B° — >■ D 7i + ir~ decays may, 
however, give contributions with distinctive shapes in the mass distributions of Dnir and 
DKtt candidates, respectively, and are therefore explicitly modelled in the fits. 

The Dtttt fit includes a double Gaussian shape to describe the signal, where the 
two Gaussian functions share a common mean, together with an exponential component 
for partially reconstructed background, and a probability density function (PDF) for 
A b — >• D°piT + decays. This PDF is modelled using a smoothed non-parametric function 
obtained from simulated data, reweighted so that the D°ir + invariant mass distribution 
matches that observed in data. The shape of the combinatorial background is essentially 
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Figure 2: Fits to the B9-, candidate invariant mass distributions for the (a) Dtttt and (b) DKtt 
samples. Data points are shown in black, the full fitted PDFs as solid blue lines and the 
components as detailed in the legends. 



linear, but is multiplied by a function that accounts for the fact that candidates with 
high invariant masses are more likely to fail the B^ mass constrained fit. There are ten 
free parameters in the Dtttt fit: the double Gaussian peak position, the widths of the two 
Gaussian shapes and the relative normalisation of the two Gaussian functions, the linear 
slope of the combinatorial background, the exponential shape parameter of the partially 
reconstructed background, and the yields of the four categories. The result of the fit to 
the Dtttt candidates is shown in Fig. [2(a) and yields 8558 ± 134 B° — > D°tt + tt~ decays. 

The DKtt fit includes a second double Gaussian component to account for the presence 
of both B° and B° s decays. The peaking background PDF for B° — > D°tt + tt~ decays is 
modelled using a smoothed non-parametric function derived from simulation, reweighted 
in the same way as described for A b — > D°pTT + decays above. The dominant partially 
reconstructed backgrounds in the DKtt fit are from B° s decays and these extend into 
the B° signal region. Instead of an exponential component, a background PDF for 
B® — > D*°K~tt + decays is included, modelled using a smoothed non-parametric func- 
tion obtained from simulation. Studies using simulated data show that this function can 
account for all resonant contributions to the B® — > D*°K~tt + final state. The function 
describing the combinatorial background has the same form as for the Dtttt fit. The 
DKtt fit has eight free parameters; the parameters of the double Gaussian functions are 
constrained to be identical for the B° and B® signals, with an offset in their mean values 
fixed to the known B°-B® mass difference [32J. The relative width of the broader to the 
narrower Gaussian component and the relative normalisation of the two Gaussian func- 
tions are constrained within their uncertainties to the values obtained in simulation. The 
result of the fit is shown in Fig. [2(b) and yields 815 ± 55 B° -> D°K + tt- and 2391 ± 81 
B® — > D°K~tt + decays. All background yields in both fits are consistent with their ex- 
pectations within uncertainties, based on measured or predicted production rates and 
branching fractions and background rejection factors determined from simulations. 



5 



4 Calculation of branching fraction ratios 



The ratios of branching fractions are obtained after applying event-by-event efficiencies 
as a function of the Dalitz plot position. The branching fraction for the B° — y D°K + -k~ 
decay is determined as 

_ B (B° -> D°K+ix~) _ N COTT (B° D°K+7r-) 
B ° = B (B° -> DH+TT-) ~ N co ^(B° D°n+ir-) ' ^ ^ 

and the branching fraction of the B® — y D°K~tt + mode is determined as 

R _ g (g z: ^° A ^ + ) _ ^ V 1 -> D°/i-7r+) 

B ° £ (5° -»• D tv+tv-) \f d J N co "(B° -> DOvr+vr-) ' U 

where the efficiency corrected yield is j\r corr = ^ PVj/e* ot . Here the index i runs over all 
candidates in the fit range, Wi is the signal weight for candidate i, determined using the 
procedure described in Ref. [28], from the fits shown in Fig. [2] and e* ot is the efficiency for 
candidate i as a function of its Dalitz plot position. The ratio of fragmentation fractions 
is f s /fd = 0.256 ± 0.020 (31] • The statistical uncertainty on the branching fraction ratio 
incorporates the effects of the shape parameters that are allowed to vary in the fit and 
the dilution due to event weighting. Most potential systematic effects cancel in the ratio. 

The PID efficiency is measured using a control sample of D*~ — > D tt~, D° — > K + ir~ 
decays to obtain background-subtracted efficiency tables for kaons and pions as a function 
of their p and pp The kinematic properties of the particles in signal decays are 

obtained from simulation in which events are uniformly distributed across the phase space, 
allowing the PID efficiency for each event to be obtained from the tables, while taking 
into account the correlation between the p and px values of the two tracks. The other 
contributions to the efficiency (detector acceptance, selection criteria and trigger effects) 
are determined from phase space simulation, and validated using data. All are found to be 
approximately constant across the Dalitz plane, apart from some modulations seen near 
the kinematic boundaries and, for the DKtt channels, a variation caused by different PID 
requirements on the pion and the kaon. The efficiency for each mode, averaged across 
the Dalitz plot, is given in Table [1] together with the contributions from geometrical 
acceptance, trigger and selection requirements and particle identification. 

The Dalitz plots obtained from the signal weights are shown in Fig. El The 
B° — > D°7t + 7t~ plot, Fig. |3]^a), shows contributions from the p°(770) and /2(1270) reso- 
nances (upper diagonal edge of the Dalitz plot) and from the D2~(24Q0) state (horizontal 
band), as expected from previous studies of this decay [T2JH3] - The B° — > D K + n~ plot, 
Fig. [3(b), shows contributions from the i^*°(892) (upper diagonal edge) and from the 
Z^2~(2460) (vertical band) resonances, also as expected [TU]. The B® — > D°K~it + plot, 
Fig. [3(c), shows contributions from the i^*°(892) (upper diagonal edge) and from the 
_D*2~(2573) (horizontal band) states. The former contribution is as expected [11]. The de- 
cay B° s — y D*2 (2573)tt + has not been observed previously but is expected to exist given 
the observation of the B® —y D* s2 ~ (2573){i + vX decay [36] . 
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Table 1: Summary of the efficiencies for Dtttt and DKtt in phase space simulation. Contri- 
butions from geometrical acceptance (e geom ), trigger and selection requirements ( e tri g &sel ) and 
particle identification (e PID ) are shown. The geometrical acceptance is evaluated for B mesons 
produced within the detector acceptance. Values given are in percent. 

B° -)> Dtttt B° -> DKtt B°, -> DKtt 
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Figure 3: Efficiency corrected Dalitz plot distributions for (a) B° — > D°7t + tt , 
(b) B° — s> D°K + tt~ and (c) B® — > D°K~tt + candidates obtained from the signal weights. 



5 Systematic uncertainties and cross-checks 

Systematic uncertainties are assigned to both branching fraction ratios due to the following 
sources (summarised in Table [2]). Note that all uncertainties are relative. The variation 
of efficiency across the Dalitz plot may not be correctly modelled in simulation. A two- 
dimensional polynomial is used to fit the variation across the Dalitz region of each of the 
four contributions to the efficiency (detector acceptance, selection criteria, PID and trig- 
ger effects). These polynomials are used to generate 1000 simulated pseudo-experiments, 
varying the fit parameters within their uncertainties. Each set of simulations is used to 
calculate the efficiency corrected yield. The standard deviation from a Gaussian fit to 
these yields is used to provide a systematic uncertainty for each decay mode. This leads 
to a systematic uncertainty of 3.4% (3.1 %) for Rbo (Rb°)- The DKtt fit model is varied 
by scaling the signal PDF width ratio to account for the different masses of the B° and B° s 
mesons, replacing the PDFs of the background components with unsmoothed versions, 
adding components for potential background from B° s — > D*°K*° and A b — > D°piT + de- 
cays, and replacing the double Gaussian signal components with double Crystal Ball [3T] 
functions. The Dtttt fit model is varied by replacing the PDF of the A b — > D°pTT + compo- 
nent with an unsmoothed version, varying the slope of the combinatorial background and 
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Table 2: Systematic uncertainties on R B o and R B o. The total is obtained from the sum in 
quadrature of all contributions. Note that all uncertainties are relative. 



Uncertainty (%) 



Source 


dO 
B v 


7^0 

B s 


Modelling of efficiency 


3.4 


3.1 


Fit model 


6.3 


4.3 


D*± veto 


< 0.1 


1.0 


D ± veto 


2.0 


0.2 


Df veto 


0.2 


0.5 


Trigger 


1.0 


1.0 


Particle identification 


1.0 


1.0 


Simulation statistics 


2.0 


2.0 


fs/ fd 




7.8 


Total 


7.8 


9.8 



replacing the exponential partially reconstructed background component with a PDF for 
B° — > D*°7i + 7T~ decays. Combined in quadrature, these contribute 6.3% (4.3%) to R B o 
(i? B o). Variations in the D^, D ± and Df vetoes contribute to R B o (Rb°), at the level of 
<0.1%, 2.0% and 0.2% (1.0%, 0.5% and 0.2%), respectively. In addition, the possible 
differences in the data to simulation ratios of trigger and PID efficiencies between the two 
channels (both 1.0%) and the limited statistics of the simulated data samples used to 
calculate efficiencies (2.0%) affect both R B o and Rb°- The uncertainty on the quantity 
fs/fd (7.8%) affects only R B o. The total systematic uncertainties are obtained as the 
quadratic sums of all contributions. 

A number of cross-checks are performed to test the stability of the results. Based upon 
the hardware trigger decision, candidates are separated into three groups: events in which 
a particle from the signal decay created a cluster with enough energy in the calorimeter 
to fire the trigger, events that were triggered independently of the signal decay and those 
events that were triggered by both the signal decay and the rest of the event. The data 
sample is divided by dipole magnet polarity. The neural network and PID requirements 
are both tightened and loosened. The PID efficiency is evaluated using the kinematic 
properties from D°ir + ir~ data instead of from simulation. The requirement for the B^ 
mass constrained fits to converge is removed. All cross-checks give consistent results. 



6 Results and conclusions 

In summary, the decay B° s — > D°K~ii + has been observed for the first time, and its 
branching fraction relative to that of the B° — >• D tt + 'k~ decay is measured to be 

B (B° D°K~7r + ) 

y = f = 1.18 ± 0.05 stat. ± 0.12 syst. . 

B (B° DH+ti-) v ' v ' 
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The current world average value of B (B° -> D°n + n-) = (8.4 ± 0.4 ± 0.8) x 10~ 4 
assumes equal production of B + B~ and B°B° at the T(4S) resonance and uses 
the D° branching fraction B (D° — >■ A'~7r + ) = (3.80 ± 0.07) %. Using the current 
world average values of T(T(45) -> £ + £-)/r(r(4S) 5°^°) = 1.055 ± 0.025 [22] and 
S (-D — > K~n + ) = (3.88 ± 0.05) % [32], the branching fraction of the normalisation chan- 
nel becomes B (B° — > D°7r + 7T~) = (8.5 ± 0.4 ± 0.8) x 10" 4 . This corrected value gives 

B (B° s D°K~n + ) = (1.00 ± 0.04 (stat.) ± 0.10 (syst.) ± 0.10 (B)) x 10~ 3 , 

where the third uncertainty arises from B (B° — > D°n + n~). The B° — > D°K + n~ decay 
has also been measured, with relative branching fraction 

B (B° -)• D°K + tv-) 

y = f = 0.106 ± 0.007 stat. ± 0.008 syst. . 

B (B° -)• DH+71-) K J v ' 

Using the corrected value of B (B° — > D°7r + ir~) gives 

B (B° D°K + n-) = (9.0 ± 0.6 (stat.) ± 0.7 (syst.) ± 0.9 (B)) x 10" 5 , 

which is the most precise measurement of this quantity to date. Future studies of the 
Dalitz plot distributions of these decays will provide insight into the dynamics of hadronic 
B decays. In addition, the B° — > DK + ir~ decay may be used to measure the CP violating 
phase 7. 
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